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Abstract 



Electron tunneling in ferromagnetic single-electron transistors is consid- 
ered theoretically in the sequential tunneling regime. A new formalism is 
developed, which operates in a two-dimensional space of states, instead of 
one-dimensional space used in the spinless case. It is shown that spin fluctu- 
ations can be significantly larger than the charge fluctuations. The influence 
of discrete energy spectrum of a small central electrode on tunneling current, 
charge and spin accumulation, charge and spin fluctuations, and on tunnel 
magnetoresistance is analyzed in details. Two different scales are found in 
the bias dependence of the basic transport characteristics; the shorter one 
originates from the discrete energy spectrum and the longer one from discrete 
charging of the central electrode. The features due to discrete spectrum and 
discrete charging disappear at high temperatures. 

73.23.Hk, 73.40.Gk, 75.70.-i 
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I. INTRODUCTION 



Electron tunneling in ferromagnetic junctions is of current interest due to expected ap- 
plications in magnetic storage technology and in other spin-electronics devices. Most of 
experimental and theoretical works published up to now deal with tunnel magnetoresis- 
tance (TMR) in simple planar junction, i.e., with variation of the junction resistance when 
magnetic configuration of the junction is changed. Tunneling in more complex junctions, 
particularly in mesoscopic ones, where charging effects are important, was studied only very 
recently. A special kind of such junctions are double junctions with a small central electrode 
(called alternatively island in the following). Tunneling in such junctions, known also as 
Single Electron Transistors (SET's), was extensively studied in the past decade, but only in 
the nonmagnetic limit. 1 It was shown that when the electrical capacitance C of the central 
electrode is small enough, the charging energy E c = e 2 /2C can be larger than the thermal 
energy k B T and discrete charging of the central electrode can lead to Coulomb blockade of 
electric current below a certain threshold voltage and to characteristic 'Coulomb staircase' at 
higher voltages. However, the interplay of ferromagnetism and discrete charging was studied 
only very recently. 2-7 It has been shown that discrete charging can lead to oscillations in 
TMR. 4 In Ref. [4] the intrinsic spin relaxation time on the central electrode was assumed to 
be sufficiently short (of the order of the time between successive tunneling events or shorter) 
to neglect spin accumulation. Apart from this, quantization of energy levels of the central 
electrode was neglected and the considerations were restricted to the limit where orthodox 
tunneling theory is applicable, i.e., to the case where the barrier resistances are larger than 
the quantum resistance Rq, Rq = e 2 /h. In that limit higher order processes (cotunneling) 
can be generally neglected, except in the Coulomb blockade regime, where they can play 
an important role and can significantly enhance TMR. 6 When the intrinsic spin relaxation 
time on the central electrode is sufficiently long (much longer than the time between suc- 
cessive tunneling events), spin accumulation on this electrode has to be taken into account 
and can lead to new phenomena. 5 First, the spin accumulation can enhance TMR. It can 
also generate TMR when the central electrode is nonmagnetic. Second, it can give rise to a 
negative differential resistance. Third, it can reverse sign of the tunnel magnetoresistance. 

Quantized nature of energy spectrum of a small central electrode and fluctuations in the 
spin accumulation were ignored in the works on magnetic SET's done up to now. These 
restrictions are relaxed in the present paper, where both energy level quantization and spin 
fluctuations are taken explicitly into account. Some preliminary results have been published 
elsewhere. 8 Accordingly, we consider a double junction in which all three electrodes can 
be ferromagnetic. In a general case relative orientation of magnetic moments of the three 
electrodes can be arbitrary. When the three electrodes have different coercive fields, shape 
anisotropy and/or some of them are exchange biased, then the magnetic configuration can be 
easily controlled by a small external magnetic field. However, we restrict our considerations 
to the case where the magnetization of one of the external electrodes and of the island 
are parallel to one another and parallel (parallel configuration) or antiparallel (antiparallel 
configuration) to the magnetization of the second external electrode. General geometry of 
the junction considered in this paper is shown schematically in Fig.l. 

In Section 2 we describe the formalism used for calculating electric current, junction 
resistance and other characteristics of the system. Numerical results are presented and 
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discussed in Section 3. Summary and final conclusions are in Section 4. 

II. DESCRIPTION OF THE METHOD 

The formalism described in this section is a generalization of the formalism developed 
for spinless SET's. 9 ' 10 ' 11 We consider a double junction in which the external electrodes are 
ferromagnetic, while the central one is either magnetic or nonmagnetic. The junction is 
shown schematically in Fig.l, where spin dependent discrete energy levels of the central 
(magnetic) electrode are also indicated. When the central electrode is nonmagnetic, the 
energy levels are spin degenerate. Generally, we assume that the left and central electrodes 
have parallel magnetizations, while the magnetic moment of the right electrode can be 
changed from antiparallel to parallel alignment (e.g. by applying an external magnetic field), 
as indicated in Fig.l. A bias voltage V is applied in such a way that the right (left) electrode 
is the source (drain) electrode for electrons. A gate voltage Vq is applied capacitively to 
the central electrode (not shown in Fig.l). Apart from this, we assume that electron spin is 
conserved during tunneling through the barriers and the spin dependent resistances of the 
left (Ri a ) and right (R ra ) junctions are larger than the quantum resistance Rq. 

Let E ia denote the single-electron energy levels of the central electrode at V = 0. The 
index if runs over all energy levels for spin a =f, while the index runs over all energy levels 
for spin u =[. The discrete energy levels include contributions from all magnetic and 
nonmagnetic interactions within the central electrode, like electron correlations responsible 
for ferromagnetism, magnetic anisotropy, etc (the Zeemann term is neglected as the magnetic 
field assumed to control magnetic configuration is assumed to be small). Generally, the 
discrete levels depend on the number of electrons in the central electrode and on their 
distribution. In our description, however, we simplify the problem and assume that the 
discrete levels are independent of the electron distribution, so the energy spectrum moves 
'rigidly' up or down when a bias voltage is applied and/or when the central electrode becomes 
charged with a certain number of excess electrons. 10 This approximation is reasonable when 
the total number of electrons on the central electrode is significantly larger than the number 
of excess electrons and larger than the number of spins accumulated on the centeral electrode. 

When a bias voltage V is applied, then a stationary electric current flowing through the 
junction is then given by 

" ia {n} 

-5K,0]/K + £j+-£ f )} = -eEEEWW) 

CT i CT {n} 

x{8[n ia , 1][1 - f{E ia + E r x* - E F )\ - 5[n ia ,0]f(E ia + E r + - E F )} , (1) 

where {n} denotes a particular distribution of the occupation numbers, {n} = {rif,n^} = 
{ni v ...,n iv ...;ni v ...,n iv ...}, of the energy levels E iiy , with n ia = 1 (n ia = 0) when the 
energy level %„ is occupied (empty). P({n}) is the stationary probability of the configuration 
{n} while 8[n, n'] is defined as 8[n, n'] — 1 for n — n' and 8[n, n'] = for n ^ n'. Apart from 
this, e denotes the electron charge (e > 0), E 1 ^* and E^z are defined as E 1 ^* = eV^* ± E c 
and E r ^z = —eV^* ± E c , where E c = e 2 /(2C) is the charging energy, N* is the number of 
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excess electrons on the central electrode, and V^* (V^+) is the electrostatic potential drop 
on the left (right) junction, 

yi __ C r + Cg N*e Co 

v n* - — q — v + —g v g , {*) 



V',, = §V - ^ + ^V . (3) 

Here, C\ and C r denote capacitance of left and right junctions, respectively, Cq is the gate 
capacitance, and C is the total capacitance of the central electrode, C = C\ + C r + Cq. 
When writing Eq.(l) we also assumed Fermi-Dirac distribution function, f(E — E F ), of the 
charge carriers in the external electrodes, with E F denoting the Fermi level (as in Fig.l). 
Finally, Y\ a (T r ia ) in Eq.(l) is the tunneling rate of electrons from the left (right) electrode 
to the level E ia of the island, 

rf) = ^|Mf } | 2 ^ {r) , (4) 

where M 1 ^ is an average matrix element for transitions from the left (right) electrode to 
the level i a and D 1 ^ is the spin dependent density of electron states in the left (right) 
electrode. We assumed above that the charging energy E c is independent of the number 
of electrons on the central electrode and on their distribution. This is usual approximation 
within the 'orthodox' description of single electron tunneling. This approximation is valid 
for thermalized distribution of electrons in the central electrode. When the electrons on 
the central electrode are not in thermal equilibrium, then the charging energy depends 
on a particular distribution of the electrons, as shown recently, both experimentally 12 and 
theoretically. 13 Taking into account the assumption of partial thermalization of electrons at 
the central electrode, as will be described later, we assume E c to be constant. 

The number of electrons with spin a on the central electrode is equal N a = J2i a n v an d 
the total number of electrons is N = N^ + Ni. It is convenient for future analysis to introduce 
also the number of excess electrons of a given spin orientation a as, N* = N a —N 0a , where N 0a 
is the number of electrons with spin a in equilibrium (at V = 0). Note that iV* = + iVj\ 
Magnetic moment of the island is then determined by the number M — — N^, while the 
excess magnetic moment by the number M* = M — M , where M is the equilibrium value 
of the number M at V = 0, M = iV 0T - iV^. 

The probability P({n}) can be determined from a stationary solution of the following 
master equation: 

+ 5^P({ni t ,...,n( i _i) T ,7i iT = l,n (i+ i) T ,...;n i })B(i T |{n}) 
+ J2 P ({ n fi n h,-, n {i-i)^ n ii = !» ■■■}) B ( i l\{ n }) 
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+ J2 p (i n i v -,n(i_i) T ,ni T = 0,7i( i+ i) T , ...;nJ)C(z T 

+ 53P({n T ;7ii i ,...,ri( i _i) i ,n il = 0, n^+i^, ...})C(i l 
h 

+ £ P({^i T , ...,ra (i _i )T ,ra iT = l,n( i+ i) T ,...,ny_i) t ,n iT = 0, ny+i) T , ...; nJ)i}(i T , j T 
+ 2 P (W; n ip = ij^i+i^, ...,n {j _ 1)l ,n jl = 0,n {j+1)l , ...})D(i l ,j l 



\{n}) 

IW) 
IW) 
IW) 
IW) 



+ [ P ({ n i T ' -' n (*-l)T' n iT = !>«(i+i)T» = °' n (i+i)i' •••}) 5 '(n ) J|IW) 



+P({ni T ,...,n ( i_i) T ,n iT =0,n( 1+1)t ,...;ni 1 ,..,n (H)l ,n Jl = 1, ny+ij^, -})^0'i, «t I W)] > ( 5 ) 
where 

A(i„|{n}) = % v , 0] {r( CT /(E jCT + - -Ef) + ^/(^ + E r + - E F )} 

+5K, i] {r*ji - f{E ia + - e p )\ + r^[i - f{E ia + ^ - # F )]} , 

B(i a \{n}) = 8[m., 0] x jrUl - /(^v + tffc +a - Sjp)] 



(6) 



+r^[i-/(^ CT + ^ +1 -^)]| 

C(v| W) = % v , 1] {^/(^ + E l +_, - E F ) + rj(E la + - £ F )} , 

#(W>IW) = ^o^K^O] 5[n iCT , 1] rf j(Tii<T + 5[cr, -<7']5[n i(T , 0] 6[nj_ a , 1] Wy_ CT ,i CT , 

D(ia,ja\{n}) = 5[n ia ,0] 5[n ja , 1] d jaiia , 

j-alW) = 5 K,>°] %i-a> 1] • 



(7) 



(9) 
(10) 

(11) 



The first term in Eq.(5) describes the rate at which a given distribution decays due to 
electron tunneling to and off the central electrode. The second and third (fourth and fifth) 
terms, on the other hand, describe the rate at which the probability of a given distribution 
increases due to tunneling processes from (to) the central electrode to (from) the external 
ones. The terms with the coefficients H, D and 5* describe the electronic (spin-conserving) 
and spin-flip relaxation processes inside the central electrode. The transition probability 
from the level Ei to the level E» is gL and to the level Ej is Wi ,■ . The master 
equation ([5]) has a general form which includes internal relaxation processes on the island 
and also the influence of gate voltage Vq- In Eq.(5) we assumed that electrons in the source 
and drain electrodes are in thermal equilibrium, while the electrons in the island can be 
generally out of equilibrium. 

It is convenient to define the probability P(N^, N±), that the island is occupied by 
electrons with spin o =f and iVj_ electrons with spin o =|, respectively. 



P(N h N t ) = £ P({n})5[N h £ ni T ] S[N it £ 



11; 



(12) 
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^From Eq.(^|) one finds the following master equation for P(Nf, iVjJ in the stationary state: 

= dP{N ^ Nl) = -P(N h N±)[A(Ni, N{) + H(N h N t )] 

+P(iV T + 1, N^B^ + 1, Nfi + P(N h N t + l)S|(iV tj Ni, + 1) 
+P(iV T - 1, A^C^ - 1, N t ) + P(tf T , - l)Q(iV T , ^ - 1) 
+P(iV T + 1, ^ - 1)S U (% iVj_) + P(iV T - 1, iV 4 + l)5 iiT (iV T , iVJ . (13) 

We have defined here the following parameters: 

A(N h Ni) = E Et 1 " F{E ia \N h N\)]{r l i<r f(Ei a + E% - E F ) + VlJ{E ia + E r + - E F )} 

+F(E i jN h iV x ){rUl - /(^ CT + JSfe - + a[l - /(£v + ^ - #*■)]} , (14) 

+r^[i-/(p v + p--p F )]}, (is) 

C a (N h Ni) = E[i - F(FjiV T , iV i )]{rt/(P lCT + Pjv + * - E F ) + T[ ct /(P 1ct + E$ -E F )}, (16) 

^^^EEE^^i^W^I^^ > ( 17 ) 
^(^^^E^^i.^^-^J^.*^ » ( 18 ) 

and also have introduced the function 

f(^|jv t ,jv 4 ) = PfJ y 1 J yN £Wk.i] 5[iv T ,E^ T ] W,E^L (is) 

V T' 1/ {n} if jj 

which is the probability that the level E ia is occupied when the island contains electrons 
of spin a =f and iV 4 electrons of spin a ={. Note, that in Eqs (14) to (16) N* is the number 
of excess electrons on the island corresponding to the numbers iVj and iVj_. When either iVj 
or iVj_ increases (decreases) by one, the corresponding number N* also increases (decreases) 
by one. 

One can easily show that the following relations are fulfilled: 

i-f(^|jv Tj jv 4 ) = pJ jy) £ f (WMK.o] <W,E^t]<W>£^-J > ( 2 °) 

V T' I) { n } i T j l 

F(E 3 _jN h N{) [1 - F(EjN h NJ] = 1 E P ({n}) 

^( iV T> iV iJ {n} 



X 



*K,0] 5[ nj _ a ,l] S[NhE^n N hJ2 n n] ■ (21) 



j i 



In the following we restrict ourselves to the case of short electronic (spin-conserving) 
relaxation time, di a ^ a 3> Tl ,T\ , while the spin relaxation time is much longer, c£j CT , jCT ^> 
Fi , T ■ ^> Wi v j r The fast electronic relaxation leads to thermalization of electrons with 
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a given spin orientation. The two spin subsystems, however, are not in equilibrium and 
correspond to different chemical potentials /ij and which are determined by TVf and 
respectively. 

The free energy of internal degrees of freedom can be expressed as 



F(N h N l ) = -k B Tj2^ 



<W, E n vl ex p ( E 



(22) 



and the probability F(E i(T \N^ iVjJ is then given by the following expression 
F^jA^iVj = exp f ^f'^ ) E^K, 1] <W,E™nl 

x^E^Jexp ^- J-^^S^j • (23) 

In the limit k B T ^> AE the distribution function F can be approximated by the Fermi- 
Dirac distribution 

F{E la \N h N{) = f(E ia - pt a (N a )) , (24) 
where the chemical potential /i a (N a ) is to be determined from the equation 

J2f(E ia -pL a (N a )) = N a . (25) 

In the regime k B T < AE, the distribution function for only two levels is significantly different 
from zero or one, so one may treat the system as effectively a two-level one. 10 If we denote 
the relevant energy levels as E 1<r and E 2a , then from the Gibbs distribution one finds the 
following expression for the function F: 

-l 



771/ ,, r , M exp(-EiJk B T) J 

exp(-Ei (T /fc B T) +exp(-E 2(r /fcBT) { 



E ia - iCgi^ + E 2a 
\k B T 



(26) 



for v = \ a and « CT = 2 a . 

When we express electric current / (see Eq.(l)) in terms of the distribution function F, 
then it is given by 

' = * E T,T,nN h N l ){[l-F(E i jN h N l )]rif(E ia +E r N t-E F ) 

N- [ ,N l (J v 

-FiEM^N^riil - f(E ia + E^ - E F )}}. (27) 

The assumption of thermal equilibrium (for a particular spin orientation) on the central 
electrode requires r in < tj, where Tj = e/I is the injection time and r in is the inelastic 
relaxation time. 10 Thus, our further analysis is valid when r in < tj. At low temperatures 
the main contribution to Tj n is due to electron-electron and electron-phonon interactions. 
Experimentally, Tj„ was extensively studied in the past by means of the weak localisation 
phenomenon and typical values of r in were found to be between 1CT 12 and 1CT 10 sec. 17-19 In 
small clusters the relaxation time r in can be larger than its corresponding bulk value. 13 
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III. CHARACTERISTICS OF FERROMAGNETIC SET'S 



In this section we describe numerical results obtained on basic characteristics of the 
junction. To simplify the picture arising from discretization of energy levels of the island, 
we assume that the levels are spin degenerate (nonmagnetic island) and equally separated 
with the inter-level spacing AE. In that case the numbers iVot and Nqi are equal, so that the 
excess magnetic moment is equal to the total magnetic moment, i.e., M* = M. Assuming 
additionally that the density of states _D^ r ) in external electrodes and the matrix elements 
M l £p are independent of energy (M 1 ^ =M l J- r ^) , one can rewrite Eq.(27) as 




-F{E ia \N h 7VJ[1 - f(E ic + E r N ~ - E F )]} , (28) 

where R r a is the resistance of the right junction, (R^)" 1 = (27r/h)\M£\ 2 D^(l/AE). Intro- 
ducing in a similar way also the resistance R l a of the left junction, one can express the 
parameters (14) to (18) in terms of R r a and R l a and then calculate the probability P(N^, N±) 
from the master equation (13). 

The formalism described above makes use of the two-dimensional space of states (jVj, N±), 
in contrast to the spinless case, where the relevant space is one-dimensional. Basic physical 
characteristics of the system are then determined by the probability P(N^,N±) introduced 
in Eq.(ll). When expressed in terms of N* and iV*, this probability will be denoted as 
P*(N*, N*). In the (N*, N*) space this probability is localized on a small number of points, 
as shown in Fig. 2 for a few values of the bias voltage in the parallel and antiparallel configu- 
rations. The area of the black dots located at the points in the (N*, N*) space is proportional 
to the corresponding probabilty P*(iV*, N*). In each case the total area of all black dots is 
normalized to unity. For V = 5 mV there is no excess electron on the island (N* = N* = 0), 
as this value of V is within the Coulomb blockade region (see Fig. 3a, where the correspond- 
ing I — V curves are shown). For V = 20 mV (within the first plateau above the threshold 
voltage in the I — V curves) the points with dominant probability P*(N*,N^) are located 
on the line corresponding to N* + N* = N* = 1. There are also points corresponding to 
N* = 0, but the corresponding probabilty is significantly smaller and these points will be 
neglected in the further discussion. It is interesting to note that different points correspond 
to different values of the excess spin on the island. In the parallel configuration these points 
are distributed symmetrically on both sides of the line corresponding to Nf = N*. Con- 
sequently, the average spin accumulated on the island is zero, contrary to the antiparallel 
configuration, where the average spin accumulated on the island is nonzero. For V = 30 mV 
the probability P*(N^,N*) is significant for N* — 1 and iV* = 2. This value of V corre- 
sponds to the transition between the first and second steps in the I — V curves. As before, 
the average spin accumulated on the island vanishes in the parallel configuration, whereas 
in the antiparallel configuration it is different from zero. Note, that the number of different 
values of the excess spin on the island is now smaller. The situation for V = 40 mV is qual- 
itatively similar to that for V — 20 mV, but the number of black dots is larger. Generally, 
one can note from Fig.2, that when the bias voltage V increases, the localization area of the 
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probability P*(N*, N*) shifts to new stationary points and embraces more and more points 
in the (N*, Nf) space. Similar tendency can be observed when the temperature increases. 
Therefore, in order to get convergence in numerical calculations, the number of states taken 
into account was dynamically changed with increasing bias voltage and temperature. 

A. Bias voltage characteristics 

Figure 3a shows the current-voltage characteristics of a junction with a nonmagnetic 
island and ferromagnetic source and sink electrodes. The single-junction resistances in the 
parallel configuration have been assumed to be = 200 MQ, Riy = 100 Mf2 for the left 
junction and R r ^ = 4 MQ, R r ^ = 2 MQ for the right one. In the antiparallel configuration 
the magnetization of the right electrode is reversed and the corresponding resistances are 
i? r | = 2 Mf2 and R r ^ = 4 MQ. Note, that the same spin asymmetry factor pj = RjjR^ has 
been assumed for both junctions in the parallel configuration, pi = p r = 1/2. Owing to a 
large difference between the resistances of the left and right junctions, the Coulomb steps in 
the I-V characteristics are clearly seen. 11,17 Since C\ > C r in the case considered here, the 
threshold voltage Vth, below which the current is blocked (N* = 0), is approximately equal 
to V^h = (C/Ci)(2E c + AE)/2e ps 10.2mV. The large steps in Fig. 3a correspond respectively 
to N* = 1, 2, and their length is V p « (C/Cj)(4E7 c + AE)/2e « 19.1 mV. There are also 
additional small steps of length V s = (C/Ci)(AE/e) ~ 2.7 mV, which result from discretness 
of the energy spectrum of the island and correspond to opening a tunneling channel with a 
new value of the excess spin on the island (new value of M) . Position of the steps is clearly 
seen in the dl/dV curves shown in Fig.3b for the antiparallel configuration. The large peaks 
correspond there to the Coulomb steps while the small ones to the steps due to discrete 
energy spectrum. This behavior is qualitatively similar to that observed experimentally in 
tunneling through small Al particles 18 or through Cqq molecules. 19 

The I-V curves in the parallel and antiparallel configurations are different (solid and 
dashed curves in Fig. 3a). Consequently, the corresponding resistances of the whole system 
are also different in both configurations; R p and R ap , respectively. This, in turn, results in 
tunnel magnetoresistance (TMR), which is described quantitatively by the ratio TMR = 
(R ap — R p )/Rp. 20 The bias dependence of TMR is shown in Fig. 3c. As one can see, TMR 
oscillates with increasing V with the period V p . The amplitude of the oscillations decreases 
with increasing voltage. In the limit V ^> E c j e the system can be treated as a set of ohmic 
resistors with the total resistance R~ l = (R^ + _R r j) _1 + (R^ + i? r |) _1 . In our case the total 
limiting resistances for the antiparallel and parallel configurations are respectively R ap = 
68.65 Mfl and R p = 68 MQ, which gives the asymptotic value of TMR equal approximately 
to 0.01. This value can be larger for systems with either larger spin asymmetry in the 
single-junction resistances, or smaller difference between the resistances of left and right 
junctions. 

For the parameters assumed in numerical calculations, the incoming electrons pass 
through the less resistive and more capacitive junction, while the outgoing electrons pass 
through the more resistive and less capacitive one. In that case electrons accumulate on the 
island when a bias voltage V is applied. Fig.4a presents the bias dependence of the charge 
accumulation. The steps in the curves show that the average charge < Q > accumulated 



9 



on the island is close to le, 2e, ... and is almost constant between the steps. Plot of the 
root mean square, rms(A r *) = [< iV* 2 > — < N* > 2 ] 1 / 2 , as a function of V is presented in 
Fig. 4b. The charge fluctuations are large at the steps, where a new charge channel becomes 
open, i.e. when N* — > N* + 1. Between the steps fluctuations are rather small. 

When the right and left junctions correspond to different spin asymmetry factors, then 
not only charge but also spin is accumulated on the island. For the junction assumed in 
Fig. 3 this happens in the antiparallel configuration. The plot of < M > as a function of V is 
shown in Fig. 5a. Indeed, there is almost no spin accumulation in the parallel configuration, 
whereas a significant spin accumulation occurs in the antiparallel configuration, which varies 
oscillatory-like with increasing V. The origin of the oscillatory behavior is described in Ref. 
[5], Here, we only note that beginning from the threshold voltage, the average < M > 
increases with increasing V up to < M >= 3, which occurs at V ~ 20 mV. At this value of 
V a new charge channel, corresponding to N* = 2, becomes open for one spin orientation, 
which reduces spin accumulation. The average < M > starts to increase again at V 30 
mV, and the second oscillation period in the spin accumulation begins. 

Figure 5b shows fluctuations of the induced magnetic moment on the island. Although 
there is almost no spin accumulation in the parallel configuration, the curve representing spin 
fluctuations in the parallel configuration is similar to that for the antiparallel one. Moreover, 
the fluctuations in M are even larger in the parallel configuration than in the antiparallel one, 
because the space of states available for fluctuations is reduced by the spin accumulation. 
Numerical analysis of the probability distribution P*(N*,N*) on the first Coulomb step 
in the parallel configuration shows two high maxima corresponding to opposite induced 
magnetic moments (which gives < M >= 0). The separation between the maxima increases 
with increasing voltage and then decreases when V exceeds 23 mV. Behavior of the spin 
fluctuation with increasing bias voltage resembles behavior of the average spin accumulated 
on the island. The fluctuations vary oscillatory-like with increasing V, with the same phase 
and period as the oscillations in < M >. It is also interesting to note, that every second 
peak of the spin accumulation and spin fluctuation in Fig.5a and Fig. 5b have a similar shape. 
This additional periodicity is due to variation of the ground state from the state with odd 
number of electrons to that with even number of electrons on the island (the (N*, N*) space 
is different along the diagonal corresponding to N* odd or N* even). 

The maximum current in our numerical results is of the order of 1 nA, which corresponds 
to the lowest value of the injection time, 17 ~ 5 x 10~ 10 sec. Thus, the numerical results 
are valid in the whole range of applied voltage when r in < 5 x 10~ 10 sec, which can be 
obeyed in real systems at T > 7], where T t is of the order of 1 K. This estimate is consistent 
with that in Ref. [13], where the inelastic relaxation time at 30 mK was estimated for a 
small Al cluster. When adapted to our value of AE, this estimate gives Tj n of the order 
of 10 -9 sec. Since Tj n decreases with increasing temperature, 7] = 1 K as the lower limit 
for validity of our numerical calculations seems to be quite reasonable. This temperature 
is low enough to observe the level quantization. The lowest temperature assumed in our 
numerical calculations is 2.3 K, which is above the lower limit T\ and also sufficiently below 
the upper limit, determined by the condition ksT AE, above which the quantization 
effects disappear. It is also worth to note, that tj can be made longer by an increase in the 
junction resistances. Thus, for realistic Tj n one can always find a range of parameters, where 
our description is valid. 
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B. Temperature dependence 



The numerical results presented above were calculated for /c B T much smaller than the 
charging energy E c and also smaller than the level spacing AE. The two energy scales were 
then clearly seen in all characteristics of the system. When the temperature increases the 
probability P*(N^, N*) spreads over larger area in the (N^ , N*) space and the peaks become 
smaller. This has a significant influence on transport properties. 

In Figs. 6a and 6b we show the current characteristics (a) and TMR (b) for different 
temperatures. The Coulomb steps in the I-V curves disappear at high temperatures, and the 
current becomes ohmic with the classical value of the resistance, R^ 1 = (R^ + R r ^)~ 1 + (Rii + 
Rrl) -1 . The corresponding value of TMR is then equal to about 0.01 and is almost voltage 
independent. The small steps in the I— V curves and TMR, which result from discreteness of 
the electronic structure of the island, disappear rather quickly with increasing temperature, 
much earlier than the Coulomb steps do. For T = 11.6 K (the thermal energy is equal to 
1 meV and is three times smaller than AE = 3 meV) most of the small steps disappear, but 
there are still well defined Coulomb steps and large oscillations due to charging effects (in 
this case E c = 6.02 meV). 

The influence of increasing temperature on the bias dependence of charge accumulation 
and charge fluctuations is shown in Fig. 7a. and 7b, respectively. The curves representing 
charge accumulation at different temperatures are similar to the corresponding I-V char- 
acteristics. The effects due to discreteness of energy levels and due to discrete charging 
gradually disappear with increasing temperature. At T = 58 K the charge accumulation 
becomes a linear function of V, as it should be in an ohmic system. Oscillations in the 
charge fluctuations are less sensitive to the temperature. As follows from Fig. 7b, they are 
periodic functions of V and the periodicity survives even at T = 58 K, where the I — V 
curves have already ohmic character. 

Spin accumulation and spin fluctuations at different temperatures are shown in Fig.8. 
The oscillations with increasing bias voltage disappear when the temperature increases, 
quite similarly as the oscillations in charge accumulation and charge fluctuations (Fig. 7). 
At T = 58 K the spin accumulation and spin fluctuations vary almost linearly with increasing 
bias. 

C. Gate voltage dependence 

Consider now transport characteristics of the system as a function of the gate voltage 
V g , which is related to the induced charge Q in = V g C g on the island. Assume a constant 
bias voltage which is above the threshold voltage and corresponds to the plateau between 
the first and second Coulomb steps, say V — 15 mV. Figures 9a and 9b show the I — V 
characteristics and TMR, respectively, as a function of the gate voltage, and calculated for 
T = 2.3 K, T = 11.6 K and T = 34.8 K. Figure 10, on the other hand, shows charge (a) 
and spin (b) accumulation on the island, calculated for the same temperatures as in Fig.9. 
Electric current, TMR and spin accumulation are periodic functions of V g , with the period 
V* = 2e/C g ~ 107 mV corresponding to AQ- m = 2e. The curve corresponding to charge 
accumulation is similar to that representing charge accumulation as a function of the bias 
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voltage (Fig.4a). Due to asymmetry of the states with odd and even numbers of electrons 
on the island, the period V g p is twice as long as in the spinless case. At low temperatures 
the difference between states with N* odd and even is clearly seen in all characteristics. At 
high temperatures, however, the difference between those two cases disappears and period 
becomes the same as in the splinless case. To understand this difference let us analyze the 
situation in more details. For V — 15 mV and V g = the average excess charge on the island 
is close to f e. The probability P*(N*, N*) has then large peaks for N* = Nt + Nf = 1. An 
increase in V g leads at low temperatures to an almost linear decrease of the current (Fig. 9a), 
while the charge accumulated on the island remains almost unchanged (it increases very 
slowly, see Fig. 10a). Origin of the decrease in electric current can be explained as follows. 
When V g increases the position of the Fermi level of the island shifts to lower energies. This 
effectively reduces the number of energy levels from which electrons can tunnel through the 
left junction. In our case this junction has much larger resistance than the right one, and 
therefore it is just the junction which determines electric current flowing through the system. 
Thus, an increase in V g results in a decrease in electric current. (Opposite behavior, when 
current increases with increasing V g is also possible for other parameters.) 

At V g ~ 38 mV the relevant states are (N*,N*) = (1,0) and (0,1). In the antiparallel 
configuration the probability for these states is P*(l, 0) = 0.76 and P*(0, 1) = 0.20, whereas 
in the parallel configuration P*(1,0) = P*(0, 1) = 0.48. V g = 38 mV is already close to 
the value at which ground state with one electron more on the island becomes energetically 
more convenient. A small increase of V g to 45 mV leads then to a large increase in the 
charge accumulation, A < Q >^ le. This also leads to a rapid increase in the electric 
current, roughly to the value it had at V g = 0. The current increases because from the 
transport point of view the system returns to the situation at V g = (without counting the 
discreteness of the energy spectrum). For V g = 45 mV there is only one relevant state, i.e. 
the state (1,1) with the probability P*(l, 1) = 0.70 in the antiparallel and P*(l, 1) = 0.75 
in the parallel configurations. The spin accumulation reaches then minimum at this point 
(see Fig. 10b). 

At Vg ~ 97 mV a new ground state is formed and the system goes over from the state 
with even number of electrons on the island to the state with odd number of electrons. Close 
to V g = 97 mV the relevant state is at (1,1), but a small increase of V g leads the system to a 
new stationary state, in which the states (2,1) and (1,2) are more important. This transition 
is different from the one at V g ~ 38 mV. This difference is clearly seen in TMR, which in 
this range of V g has a large deep and becomes negative. 

At higher temperatures the difference between the situations with odd and even numbers 
of excess electrons on the island is not visible (see the curves corresponding to T = 11.6 
K and 34. 8K in Fig.9 and Fig. 10). The thermal energy is then comparable to the energy 
needed for opening a new spin channel, i.e. ksT m AE. Two states in the (N* , iV*)-space, 
for which AM = ±2, are difficult to be distinguished. Therefore, the periodicity is then as 
in the spinless case. 
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IV. SUMMARY AND CONCLUSIONS 



We have developed formalism for calculating electric current, spin and charge accumula- 
tion and TMR in ferromagnetic SET's with a small central electrode - small enough so that 
the discrete structure of its energy spectrum plays a significant role. We found two different 
scales in all characteristics of the junction; the shorter one related to the discreteness of 
energy spectrum and the longer one related to discrete charging of the island with single 
electrons. The features due to discrete energy levels can be seen at low temperatures and 
disappear relatively quickly with increasing temperature; much faster than the features due 
to discrete charging (AE « E c in our case). 

The junction characteristics are periodic functions of the bias and gate voltages. At low 
temperatures the periods are twice as long as the corresponding ones at high temperature. 
This is because at low temperatures the situations with even and odd numbers of electrons 
on the island can be distinguished, while at high temperatures this difference disappears. 

We have also shown that spin fluctuations can be significantly larger than the charge 
fluctuations. Such large spin fluctuations can play a significant role in the current noise. 21 
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FIGURES 




FIG. 1. Geometry of the junction and schematic profile of the potential energy when a bias 
voltage V is applied. Discrete energy levels of the island for both spin orientations are indicated 
by the solid and dashed lines. 
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FIG. 2. Probability P*(N^, N*) in the (N^, N*) space of states (proportional to the area of the 
black dots) calculated for the antiparallel and parallel configurations and for four different values 
of the bias voltage. The junction resistances in the parallel configuration are: = 200 
Ril = 100 Mf2, i? r j = 2 MJ7 and R r ^ = 4 whereas in the antiparallel configuration 

R r ^ = 4 MJ7, R r i = 2 MQ. The other parameters assumed in numerical calculations are: C\ = 9 aF, 
C r = 1.3 aF, C fl = 3 aF, £ c = 6.02 meV, A£ = 1.8 meV, V = 26 mV, = 0, T = 2.3 K and 
^ = o. 
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V [mV] 

FIG. 3. Voltage dependence of the tunnel current / (a), derivative dl/dV (b), and tunnel 
magnetoresistance (c) determined at T = 2.3 K. The solid and dashed curves in (a) correspond 
to the antiparallel and parallel configurations, respectively, whereas the plot in (b) is for the 
antiparallel configuration only. The parameters of the system are the same as in Fig. 2. 
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FIG. 4. Charge accumulation < N* > (a) and charge fluctuations [< iV* 2 - < N* > 2 ] 1 / 2 (b) 
as a function of the bias voltage. The solid and dashed curves corresponds to the antiparallel and 
parallel configuration, respectively. The parameters are the same as in Fig.2. 
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FIG. 5. Spin accumulation < M >=< iVt — > (a) and the spin fluctuations 
[< M 2 — < M > 2 ] 1 / 2 (b) as a function of V in the system defined in Fig.2. Solid and dashed 
curves are for the antiparallel and parallel configuration, respectively. 
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FIG. 6. Bias dependence of electric current in the antiparallel configuration (a) and TMR (b) 
for different temperatures and for AE/kB = 34.8 K and E c /kB = 69.9 K. The other parameters 
are the same as in Fig. 2. 
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FIG. 7. Bias dependence of the charge accumulation (a) and charge fluctuations for different 
temperatures. The parameters are the same as in Fig. 2. 



22 




10 20 30 40 50 60 70 80 

V[mV] 



FIG. 8. Voltage dependence of the spin accumulation (a) and spin fluctuations for different 
temperatures. The parameters are the same as in Fig. 2. 
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FIG. 9. The tunneling current / in the antiparallel configuration (a) and TMR (b) as a function 
of the gate voltage V g calculated for V = 15 mV and for T = 2.3 K (solid curve), T = 11.6 K 
(dashed curve), T = 34.8 K (dotted curve). The other parameters of the system are as in Fig.2. 
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FIG. 10. The charge (a) and spin accumulation (b) as a function of the gate voltage V g for 
V = 15 mV and for T = 2.3 K (solid curve), T = 11.6 K (dashed curve), T = 34.8 K (dotted 
curve). The other parameters of the system are as in Fig. 2. 
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